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基于 PML-V2 数 据 集 的 中 国 燕 散 发 时 空 动态 特征 分 析 
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(中 国 地 质 大 学 (武汉 ) 环 境 学 院 , 湖 北 武汉 430074) 


摘 m. 陆地 藻 散 发 是 陆 - 气 水 热 交换 的 重要 过 程 ,是 全 球 水 分 运 移 及 能 量 转换 中 重要 的 一 环 。 识 
别 蒸 散发 的 时 空 动态 特征 对 研究 区 域 水 循环 和 能 量 转换 具有 重要 意义 。 此 次 研究 利用 9 个 通 量 站 
观测 数据 ,对 比分 析 了 3 种 全 球 蒸 散 发 产品 (GLEAM MOD16 以 及 PML-V2) 在 中 国 区 域 的 适用 性 。 
在 此 基础 上 , 许 取 在 中 国 区 域 适宜 性 最 佳 的 地 表 蒸 散发 产品 数据 ,对 2003 一 2020 年 中 国 地 表 蒸 散 
发 时 空 动态 特征 进行 了 分 析 。 结 果 表 明 :(1) PML-V2 蒸 散发 产品 在 中 国 区 域 的 适用 性 最 佳 。(2) 
在 时 间 上 ,中 国 蒸 散 发 在 研究 时 间 段 呈 增 长 趋势 ;空间 上 ,研究 时 间 段 由 西北 到 东南 总 体 呈 增加 趋 
势 ,(3) 基于 Hurst 指 数 分 析 的 未 来 中 国营 散发 变化 状况 ,预期 将 与 过 去 表现 相反 趋势 , 即 过 去 总 体 


Hu 


为 区 域 水 资源 利用 及 优化 配置 提供 参考 借鉴 。 
X 键 词 : ARK; 适宜 性 ; 动态 特征 ; HE 
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蒸 散 发 的 研究 对 于 水 资源 的 合理 配置 .农作物 
的 增产 增 量 以 及 区 域 气候 的 动态 分 析 意义 重大 …。 
目前 对 于 给 散发 的 直接 测量 一 般 是 集中 在 小 尺度 
范围 ,测量 方法 有 蒸 渗 仪 . 波 文 比 、 涡 度 相关 等 ,对 
于 大 区 域 尺 度 的 蒸 散发 一 般 无 法 通过 实际 测 得 ”。 
随 着 遥感 技术 和 经 验 理论 模型 的 发 展 ,使 得 大 区 域 
斥 度 的 实际 蒸 散 发 估算 成 为 了 可 能 。 目 前 全 球 尺 
度 与 区 域 拟 合 效果 较 好 的 产品 主要 有 基于 PM (Pen- 
man-Monteith ) 3X 8 MOD 16" ,PML-V2 ^ fild&-T- PT 
( Priestley-Taylor) 3X8] GLEAM? ZS BOR P" mi o 

蒸 散 发 的 物理 机 制 受 多 种 气候 环境 等 地 表 参 
量 的 影响 ,通过 对 蒸 散 发 在 全 国 时 空 动态 特征 的 分 
析 , 能 够 从 宏观 角度 对 全 国 地 域 气候 燕 散 情况 和 未 
来 蒸 散发 的 变化 情况 作 直 观 的 了 解 ,为 干旱 预警 、 
水 资源 调度 提供 数据 来 源 和 分 析 依据 。2011 年 刘 
昌明 等 5 利用 优化 后 的 PM 公 式 计 算 了 中 国 的 潜在 
蔡 散 发 ,并 分 析 了 相关 气象 环境 因子 对 菩 散 发 变化 
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增长 趋势 ,未 来 燕 散 发 会 下 降 。 此 次 研究 通过 对 中 国 地 表 燕 散发 时 空 动态 特征 进行 分 析 , 以 期 


趋势 的 影响 ; 张 巧 风 等 "基于 MODIS 蒸 散发 数据 对 
锡林郭勒 草原 近 14 a 的 蒸 散发 时 空 动态 作 分 析 ,发 
现 锡林郭勒 草原 蒸 散 发 总 体 表 现 为 由 东北 部 和 东 
南 延 线 向 西南 部 递减 旦 空间 变异 比较 稳定 ; 邓 兴 焰 
等 “利用 MODIS 蒸 散发 数据 集 对 天 山 山 区 2000 年 
到 2014 年 葵 散 发 作 了 时 空 动 态 特征 分 析 , Iz RAE BK 
发 在 空间 上 有 西部 大 东部 小 .北部 大 南部 小 , 且 近 
15 a 全 区 落 散 发 变异 程度 不 明显 ; 郑 德 凤 等 J 利用 
河北 省 12 个 气象 站 的 实际 气象 观测 资料 ,采用 PM 
公式 计算 了 河北 省 1967 一 2019 年 的 潜在 蒸 散 量 , 分 
析 发 现 研 究 时 段 内 莹 散发 存在 明显 的 地 带 性 特征 ， 
其 总 体 呈 上 升 趋势 。 对 于 蒸 散发 的 时 空 动态 分 析 ， 
大 多 数学 者 主要 集中 在 区 域 性 尺度 的 蒸 散发 分 析 ， 
或 通过 改进 PM 公式 计算 的 潜在 蒸 散发 进行 分 析 ， 
或 选取 蒸 散发 产品 分 析 时 没有 考虑 到 蒸 散发 产品 
在 区 域 的 适宜 性 。 本 文通 过 全 球 通 量 站 网 发 布 的 
FLUXNET2015 数据 集中 中 国 的 9 个 通 量 站 数据 ” ， 
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水 资源 规划 及 生态 水 文 研究 提供 参考 及 借鉴 意义 。 


1 研究 区 概况 


中 国 地 处 亚洲 东南 部 ,总 体 地 势 西 高 东 低 , 气 
候 类 型 多 种 多 样 。 根 据 中 国 科 学 院 地 理科 学 与 资 
源 研究 所 资源 环境 科学 与 数据 中 心气 修 区 划 数 据 ， 
以 日 平均 气温 不 低 于 10 稳定 期 的 积温 和 最 冷 月 
气温 或 极端 最 低 气 温 多 年 平均 值 为 热量 指标 ,以 干 
燥 度 ( 见 气候 指数 ) 为 水 分 指标 将 研究 区 划分 为 10 
个 一 级 气候 区 进行 分 区 讨论 "(图 1) ,其 中 琼 南西 
沙 区 和 南沙 区 各 岛屿 分 别 属于 中 热带 (VIID 和 南 热 
we (IX) ,岛屿 面积 较 小 且 分 散 ,对 于 总 体 研 究 意 义 


+ 通 量 站 - 8364 ， 
一 气候 分 区 B is dee! 
审 图 号 : GS(2023)549 l dd d 


图 1 研究 区 示意 图 
Fig. 1 Sketch map of the study area 


不 大 , 故 本 次 研究 不 对 2 个 区 域 的 莹 散发 作 探讨 。 
2 数据 与 方法 


2.1 数据 来 源 

2.1.1 通 量 站 数据 由 于 FLUXNET2015 (https:/ 
fluxnet.org ) 数 据 集 的 完善 性 与 缺失 值 或 异常 值 处 理 
方法 的 先进 性 ,本 次 研究 选取 数据 集中 位 于 中 国 区 
域 的 9 个 共享 数据 站 点 (剔除 了 一 个 误差 较 大 的 数 
据 站 点 ) ,对 全 球 蒸 散发 产品 在 中 国 地 区 的 适宜 性 
进行 研究 ,站 点 分 别 是 CN-Cha, CN-Cng, CN-Dan , 
CN- Din, CN- Du2, CN- Du3 , CN- Ha2, CN- HaM , CN- 
Qia ,站 点 分 布 如 图 1 ,站 点 详细 信息 如 表 1。 

2.1.2 GLEAM 蒜 散 发 产品 数据 GLEAM (Global 
land evaporation amsterdam model) Z& fI Az P^ fin AE 
集 产 品 通过 卫星 遥感 数据 估算 0.25° 分 辨 率 每 天 的 
实际 地 表 蒸 散发 和 植物 根系 土壤 水 含量 。 该 数据 
集 算法 由 Miralles 于 2011 年 首次 发 布 "中 ,2014 年 改 
进 为 GLEAM (v2)5 , 2017 修订 发 展 为 最 新 的 
GLEAM(v3) 系 列 产品 。GLEAM(v3) 在 保持 水 分 通 
量 与 之 前 版 本 大 致 一 致 的 情况 下 改进 了 与 91 个 通 
量 站 点 和 2325 个 土壤 水 分 传感器 各 层 土壤 水 分 含 
量 的 相关 系数 。 数 据 下 载 网 址 www.GLEAM.eu。 为 
了 对 产品 进行 中 国 区 域 适宜 性 分 析 , 本 次 研究 根据 
通 量 站 数据 及 产品 数据 时 间 范 围 下 载 了 2002 年 1 
月 一 2010 年 12 月 GLEAM(v3) 每 月 的 燕 散 发 数据 。 
2.1.3 MOD16 数 据 MOD16 数 据 集 是 由 美国 国家 
航空 局 通过 Aqua Fil Terra 卫星 平台 数据 通过 
MOD16 算 法 生产 的 全 球 地 表 蒸 散发 产品 数据 集 ,该 
数据 集 产 品 包含 空间 尺度 500 m 和 时 间 尺 度 8 d 的 
地 表 蒸 散发 ,覆盖 面积 为 全 球 1.09x 10* km 的 陆地 ， 


表 1 通 量 站 点 基本 信息 


Tab.1 Basic information of flux site 


站 点 ID 站 点 名 经 度 /(9) 纬度 /(°) 土地 覆盖 类 型 年 份 

CN-Cha 长 白山 128.0958 42.4025 混 生 林 2003 一 2005 
CN-Cng 长 岭 123.5092 44.5934 草地 2007 一 2010 
CN-Dan 当 雄 91.0664 30.4978 草地 2004 一 2005 
CN-Din 易 湖 山 112.5361 23.1733 Aag n E 2003—2005 
CN-Du2 多 伦 草地 116.2836 42.0467 草地 2006 一 2008 
CN-Du3 多 伦 退 化 草 多 116.2809 42.0551 草地 2009 一 2010 
CN-Ha2 海北 灌 从 101.3269 37.6086 湿地 2003 一 2005 
CN-HaM 西藏 海北 高 山 站 101.1800 37.3700 草地 2002 一 2004 
CN-Qia 千 烟 洲 115.0581 26.7414 常 绿 针 叶 林 2003 一 2005 


202307.00197v1 


chinaXiv 


涵盖 了 各 种 植被 类 型 ,该 数据 集 的 有 效 数 据 日 期 为 
2000 年 1 月 至 今 史 。MOD16 算 法 基于 PM 公 式 计 算 
的 燕 散 发 主要 分 3 部 分 :土壤 蒸 散 发 .湿润 的 冠 层 落 
散发 和 植物 干 冠 层 表面 的 蒸腾 ,通过 计算 各 部 分 莱 
散 量 的 基础 上 分 别 计算 白天 和 夜间 蒸 散发 并 整合 
到 天 尺度 。 该 产品 数据 通过 AmeriFlux 的 46 个 通 量 
站 点 数据 与 产品 数据 集 对 比 显示 ,改进 的 算法 有 效 
地 将 平均 绝对 误差 (MAE) 由 0.39 mm.d :提高 到 
0.33 mm*d"'。 为 了 对 产品 进行 中 国 区 域 适宜 性 分 
析 , 本 次 研究 根据 通 量 站 数据 及 产品 数据 时 间 范 围 
下 载 了 2002 年 1 月 一 2010 年 12 月 每 8d 的 蒸 散 发 数 
据 ,并 通过 取 月 均值 的 方法 ,将 时 间 尺 度 处 理 为 每 
月 。 数 据 通过 GEE(Google Earth Engine ) 平 台 下 载 。 
2.1.4 PML-V2 数 据 PML-V2 数 据 集 是 中 国 科 学 院 
地 理科 学 与 资源 研究 所 Zhang 等 “12019 年 利用 谷歌 
地 球 引 擎 GEE 平台 ,以 遥感 数据 为 驱动 ,基于 Pen- 
man-Monteith-Leuning 蒸 散发 模型 计算 的 全 球 地 表 
蒸 散发 数据 ,该 数据 集 的 空间 分 辩 率 为 300 m ,时 间 
分 辨 率 为 8 d5。 该 数据 集 通 过 全 球 91 个 广泛 分 布 
的 通 量 站 点 (涵盖 了 10 种 植被 类 型 ) 数 据 验证 ,产品 
数据 集 的 蒸 散发 均 方 根 差 和 偏差 有 效 的 减 小 到 
0.69 mm…d-!' 和 -1.8%。PMIL-V2 数 据 集 时 间 范 围 为 
2002 年 7 月 一 2020 年 12 月 ,本 次 选取 该 时 间 段 的 全 
部 数据 ,并 通过 取 月 均值 的 方法 ,获取 2002 年 7 月 一 
2020 年 12 月 的 月 尺度 数据 ,用 以 适应 性 分 析 对 比 。 
2.1.5 其 他 参考 数据 中 国 区 域 的 DEM 高 程 数据 为 
美国 国家 航空 航天 局 喷气 推进 实验 室 SRTM v3 产 
品 55 ,分 辩 率 为 30m ,通过 GEE 平 台 下 载 。 

2.2 研究 方法 

2.21 蒜 散 发 产品 区 域 适 宜 性 评价 为 了 分 析 
GLEAM , MODI6 和 PML-V2 3 种 蒸 散发 产品 在 中 
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区 域 的 适宜 性 ,本 人 研究 通过 FLUXNET2015 数据 集 
中 中 国 区 域 通 量 站 点 的 通 量 站 蒸 散 发 数据 与 3 种 蒸 
散发 产品 的 站 点 薰 散发 数据 作对 比 ,利用 数据 拟 合 
的 皮尔 逊 相关 系数 . 均 方 根 差 分 析 3 种 蒸 散 发 产品 
在 中 国 的 适宜 性 "“"。 

(1) 皮尔 逊 相关 系数 。 皮 和 尔 逊 相关 系数 广泛 用 
于 度量 2 个 变量 之 间 的 相关 程度 ,本 研究 通过 对 比 
站 点 实测 数据 与 蒸 散发 产品 数据 的 相关 性 强度 作 
为 蒸 散发 产品 在 中 国 区 域 适宜 性 的 评价 指标 ,计算 
公式 如 下 : 


n 


Di -D0 —D) 
r= i=l (1) 
de — xP Xe: -7 

i=1 i=1 
式 中 :r 为 皮尔 逊 相关 系数 ;wy 分 别 为 样本 值 ; 元 、 
了 分 别 为 两 组 数据 的 均值 ;n 为 样本 总 数 。 

(2) 均 方 根 差 。 均 方 根 差 (RMSE ) 作 为 一 种 常 
用 测量 数值 之 间 差异 的 量度 ,本 研究 通过 站 点 观测 
值 与 蒸 散发 产品 模型 值 的 RMSE 来 评价 蒸 散 发 产品 
模拟 中 国 区 域 蒙 散发 实际 值 , 计 算 公 式 如 下 ; 


n 


2 
Y obs,i B 3 s) 


RMSE= | (2) 
n 
式 中 : x 为 第 i 个 观测 数据 ; 2,9, 为 第 ;个 模拟 
数据 。 


2.2.2 时 室 动 态 将 征 分 析 燕 散发 时 空 动 态 特 征 分 
析 主 要 是 利用 各 种 分 析 方 法 ,对 蔡 散 发 在 时 间 和 空 
间 上 的 动态 变化 进行 分 析 及 预测 ,发现 其 中 的 规 
律 ,为 更 好 规划 和 进行 水 资源 配置 提供 参考 依 
je 9. ^ gp HS pel, 


R2 分 析 方 法 
Tab.2 Analytical method 


方法 


作用 


数据 间 的 离散 程度 。 


统计 变量 各 组 数据 的 变异 程度 ,统计 分 析 莹 散发 在 时 间 序 列 上 的 稳定 性 ,变异 系数 的 大 小 反应 了 像 元 各 年 份 


Theil-Sen Median 趋势 分 析 


稳健 的 非 参 数 统计 趋势 计算 方法 ,趋势 值 (8) 大 于 0 时 ,表示 正 增长 趋势 ,小 于 0 则 表示 下 降 趋 势 。 


Mann-Kendall 检验 非 参数 的 时 间 序 列 趋势 性 检验 方法 。 在 给 定 显 著 性 水 平 下 , 当 M-K 检验 值 (2) 的 绝对 值 小 于 临界 值 时 , 即 接 


赫 斯 特 指数 


受 原 假设 ,趋势 不 显著 ; 当 Z 的 绝对 值 大 于 临界 值 时 , 则 拒绝 原 假设 ,趋势 显著 。 
用 来 衡量 时 间 序列 是 否 有 长 期 记忆 的 一 个 指标 , 赫 斯 特 指数 (HH) 位 于 0 和 1 之 间 , 当 且 为 0.5 时 ,序列 是 随机 
的 , 即 没 有 自 相 关 性 ; 当 石 大 于 0.5 时 ,序列 有 很 强 的 正 相 关 性 ,未 来 趋势 与 过 去 趋势 一 致 ; 当 且 小 于 0.5 时 , 序 


列 有 很 强 的 负 相 关 性 ,未 来 趋势 与 过 去 趋势 相反 。 
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3 结果 与 分 析 


3.4 AMA MAE 
AS XC 9 MAE e ZUR EU Tir OL UA 2 , FÉ] 3 


汇总 计算 全 国 多 年 平均 蒸 散 发 。 将 多 年 平均 蒸 散 
发 分 为 3 级 (图 4 图 5): 低 蒸 散发 区 (ET<200 mm-a”), 
rP 2 BU DX (200 mm * a" &ETE400 mm a^) i ZEB 
发 区 (ET>400 mm-a)", AeA 45 E821 h P dE 


所 示 : 


向 南 呈 逐渐 增加 趋势 ,区 域内 最 大 值 一般 出 现在 水 
域 附 近 。 由 变异 系数 可 知 , 北 温带 区 (了 ) .中 亚热带 


由 图 2 可 知 ,除了 CN-Du3 站 点 外 ,各 通 量 站 点 
数据 与 CLEAM MOD16 和 PML-V2 蒸 散发 产品 数据 
都 有 较 好 相关 性 , 蒸 散发 随时 间 变 化 的 趋势 大 体 一 
致 。 由 图 3 的 Taylor 图 可 以 观察 相关 系数 .标准 差 
和 均 方 根 差 3 个 指标 的 综合 效果 ,CLEAM 5j PML- 
V2 与 观测 值 距离 较 近 ,PML-V2 除 CN-Du3 站 点 以 
外 ,3 个 指标 均 不 如 CLEAM 和 MOD16 产 品 。 综 上 ， 
PMIL-V2 各 项 评价 指标 相对 其 他 两 种 产品 均 表现 出 
一 定 优势 ,在 时 间 和 空间 尺度 上 与 站 点 蒸 散发 数据 
的 拟 合 效果 较 好 ,趋势 一 致 ,在 中 国 区 域 有 较 好 的 
适宜 性 。 因 此 ,对 中 国 区 域 的 薰 散 发 时 空 动态 特征 
分 析 选 取 了 PML-V2 产品 数据 集 。 

3.2 2003 一 2020 年 中 国 蒸 散发 空间 分 布 特征 
将 PML-V2 蒸 散发 产品 8 dx500 m 分 辩 率 数据 


EC CV) 南亚 热带 区 (VD) . 北 热带 区 (VID) .中 热带 区 
(CVIT) 空 间 蒸 散发 变异 程度 较 低 , 葵 散发 区 域 变化 
不 大 。 中 温带 区 (1) 、 南 温带 区 (HI) 、 北 亚热带 区 
CIV) 、 高 原 气候 区 (H) 变 异 程度 较 大 , 蔡 散 发 在 空间 
分 布 相差 较 大 ( 表 3)。 贫 盗 地 区 [至 少 60% 的 面积 》 
植被 不 足 10% 的 无 植被 贫 瘫 ( 沙 、 岩 、 土 ) 地 区 ] 主 要 
分 布 在 中 国 西北 地 区 ,受气 候 和 和 覆盖 类 型 的 共同 作 
用 蒜 散 发 较 低 。 在 同一 纬度 地 区 ,不 同 土地 覆 羡 
类 型 和 所 处 的 地 理气 候 类 型 区 蒸 散发 也 有 不 同 程 
度 的 差异 所。 气候 分 区 图 与 蒸 散发 多 年 均值 有 较 
好 的 一 致 性 趋势 。 各 气候 分 区 内 薰 散 发 的 最 大 值 
随 分 区 类 型 由 西北 像 东南 逐渐 增加 。 
3.3 2003 一 2020 年 中 国 蒸 散发 时 序 变 异 特征 

由 2003 一 2020 年 中 国 燕 散发 时 序 变异 系数 空 


——obs ——GLEAM ~—MOD16 ——PML-V2 
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注 :obs 为 通 量 站 观测 数据 ;CLEAM .MOD16 .PML-V2 分 别 为 3 种 产品 数据 。 下 同 。 


图 2 站 点 拟 合 图 
Fig.2 Diagram of site fitting 
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(c) CN-Dan 


o MODI6 x PML-V2 


注 : 图 中 黑 点 观测 点 即 为 obs, 其 他 点 位 与 观测 点 距离 越 近 即 产品 效果 越 好 。 


图 3 站 点 拟 合 Taylor 图 
Fig. 3 Taylor diagram of site fitting 


间 分 布 (图 6) 和 变异 系数 数据 统计 ( 表 4) 可 以 看 出 ， 
北 温带 区 (1) .中 温带 区 (ID) 、 南 温带 区 (II) 高原 气 
候 区 (H)3 个 区 变异 系数 分 别 为 0.37、0.35、0.38、 
0.34 ,在 时 间 序 列 上 的 变异 相对 其 他 区 较 大 。 与 上 
部 分 空间 分 析 的 变异 性 类 似 , 蒸 散发 时 空 变 异 较 大 
区 域 主要 为 中 国 西北 部 ,在 沿海 地 区 城市 与 较 高 的 
时 间 变 异 系数 相关 性 较 强 。 变 异 系数 较 大 区 域 大 


多 集中 在 盆地 边缘 山脉 附近 ,滨海 城市 亦 有 较 大 的 

2003 一 2020 年 中 国 全 区 蒸 散发 时 间 上 的 变化 
特征 如 图 7 所 示 。 并 对 各 区 每 年 平均 蒸 散发 进行 一 
元 线性 回归 , 拟 合 公 式 为 : ET=B, +B, xX ,其 中 
为 年 份 , B, 为 拟 合 回归 常数 项 , B 为 拟 合 回归 一 次 
项 系数 ,各 拟 合 回归 系数 如 表 5。 
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图 4 2003 一 2020 年 中 国平 均 蒸 散发 分 级 
Fig.4 Average evapotranspiration classification of China 
from 2003 to 2020 


MES 7 n] LAG H , 2003—2020 EK DX ZEE 
FEN. KAS MAK Ne Bl EXC [a] AY 
析 趋 势 一 致 ,由 西北 向 东南 区 逐渐 增 大 。 

全 区 薰 散 发 2003 一 2020 年 波动 微弱 ,为 增 大 趋 
势 ,变化 率 为 0.89 mm'.a-'。 为 了 研究 蒸 散发 在 时 间 
上 的 变化 趋势 及 趋势 程度 的 显著 性 和 可 信和 度 ,通过 
Theil-Sen median(T-S) 和 Mann-Kendall(M-K) 分 析 方 
法 ,进一步 计算 了 2003 一 2020 年 的 趋势 斜率 (T-S 分 
析 系 数 如 图 8) 和 归 一 化 检验 统计 量 (M-K 分 析 2Z 值 
如 图 9)。 

对 中 国 各 区 进行 区 域 统计 分 析 如 表 6: 

Theil-Sen median 趋势 分 析 与 一 元 线性 回归 分 
析 基 本 一 致 ,全 区 蔡 散 发 程 增长 趋势 。 由 M-K 分 析 


人 


0 500km 
| | 

ra" " ZEB /mm:d? 

上 通 量 3564.73 jy 
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图 5 2003 一 2020 年 中 国 多 年 平均 蒸 散 发 空间 分 布 
Fig. 5 Spatial distribution of annual mean evapotranspiration 
in China from 2003 to 2020 


= 


Z 值 图 和 相应 的 区 域 统计 分 析 表 可 知 ,除了 北 温带 
区 (D) 的 置信 度 水 平 较 高 外 ,其 他 区 域 的 置信 和 度 水 平 
BEG , 即 蒸 散发 增加 或 减少 在 全 区 的 趋势 不 显著 。 
M-K 趋 势 显著 性 分 析 时 一 般 采 取 双 边 趋势 检验 , 通 
常 将 置信 和 度 90% 、95% 和 99% 对 应 Z 绝 对 值 1.65、 
1.96 和 2.58 作为 临界 值 ,将 显著 性 划分 为 : 极 显著 
(1Z1»2.58) .显著 (1.96<12I<2.58) 、 微 显著 (1.65<12ZI 
1.96) .不 显著 (0<IZI<1.65) 和 无 变化 (Z=0)。 在 中 
国 区 域内 ,2003 一 2020 年 蒸 散发 的 增长 或 下 降 一 般 
趋势 不 显著 ,但 是 在 局 部 区 域内 , 2E BUR TEE BH 
的 增长 或 下 降 趋 势 。 
3.4 ZEB B TRU 

通过 中 国 区 域 2003 一 2020 年 PML-V2 HAE Bc Az 


表 3 2003 一 2020 年 中 国 蒸 散发 统计 分 析 


Tab. 3 Statistical analysis of evapotranspiration in China from 2003 to 2020 


分 区 面积 /10' km? 最 小 值 /mm'ar 最 大 值 /mm'an 均值 /mm'.ar 标准 差 /mm +a 变异 系数 
北 温带 区 (1) 18.73 38.02 716.72 107.47 24.50 0.23 
中 温带 区 (II) 303.49 0.20 2205.41 180.85 108.76 0.60 
南 温带 区 (II1) 162.77 0.10 2466.08 259.15 215.44 0.83 
北 亚热带 区 (IV) 56.76 0.20 2283.38 537.27 232.99 0.43 
中 亚热带 区 (V) 120.64 0.20 2462.84 522.95 147.82 0.28 
南亚 热带 区 (VD) 38.23 0.10 3422.79 719.66 179.79 0.25 
北 热带 区 (VID) 7.60 0.20 3434.04 858.29 203.59 0.24 
中 热带 区 (VIII) 0.91 1.01 3382.74 990.66 256.40 0.26 
高 原 气候 区 (H) 251.41 0.10 2288.14 278.22 207.24 0.74 
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图 6 2003 一 2020 年 中 国 蒸 散发 时 序 变 异 系数 空间 分 布 


Fig.6 Spatial distribution of temporal variation coefficients 


of evapotranspiration in China from 2003 to 2020 


数据 计算 中 国 区 域 的 赫 斯 特 指数 (图 10) ,对 中 国 区 
域 赫 斯 特 指数 进行 统计 分 析 有 : 北 温带 区 (了 ) 赫 斯 特 

首 数 略微 大 于 0.5 , 即 区 域 未 来 菊 散 发 与 过 去 趋势 一 
致 , 呈 略 微 增长 趋势 。 南 温带 区 (TUI) 赫 斯 特 指数 为 
0.44 ,未 来 蒸 散发 与 过 去 呈 相 反 趋 势 , 即 略微 增长 。 
北 亚热带 区 (IV)、 中 亚热带 区 (V) 、 北 热带 区 (VID)、 
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RA 2003 一 2020 年 中 国 蒸 散 发 时 序 变 异 系数 统计 分 析 
Tab.4 Statistical analysis of temporal variation 
coefficients of evapotranspiration 
in China from 2003 to 2020 


分 区 面积 10'km 最 小 值 ” 最 大 值 ”均值 puce 
北 温带 区 (D) 18.73 0.09 128 037 0.16 
中 温带 区 (ID 303.49 0.03 424 035 0.15 
南 温带 区 (IID) 162.77 0.03 424 038 035 
北 亚热带 区 (IV) 56.76 0.00 | 424 013 0.06 
中 亚热带 区 (V) 120.64 0.00 424 0.16 0.08 
南亚 热带 区 (VD 38.23 000 424 0.13 O11 
北 热带 区 (VID) 7.60 0.03 424 0.11 0.07 
中 热带 区 (VIID 0.91 0.00» | 424 009 021 
高 原 气候 区 (H) 251.41 0.03 424 034 0.22 


中 热带 区 (VII) 、 高 原 气候 区 (H) 赫 斯 特 指数 均 小 于 
0.5 ,未 来 蒸 散 发 与 过 去 旦 相反 趋势 ,预测 未 来 蒸 散 
发 可 能 会 减 小 。 中 温带 区 (ID) 和 中 亚热带 区 (V) 由 
于 两 种 方法 趋势 性 分 析 结 果 不 一 致 ,趋势 变化 率 较 
小 , 且 置 信和 度 水 平 较 低 ,未 来 趋势 可 能 在 稳定 水 平 
波动 。 中 国 总 体 赫 斯 特 指数 小 于 0.5, 即 未 来 趋势 与 
过 去 相反 。 由 上 文 , 蒸 散发 的 时 间 变 化 特征 可 知 中 
国 蒸 散发 为 增长 趋势 ,未 来 蒸 散 发 即 可 能 与 过 去 增 
长 相反 ,逐渐 减 小 。 


-se- 北 温带 区 ( --- 南 亚热带 区 (VT) — 北 温带 区 (D) 一 南亚 热带 区 (VD 
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1200- +- 南 温 带 区 GID -4- 中 热带 区 (VIID 一 南 温带 区 (ID 一 中 热带 区 (VIID 
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注 :虚线 为 实际 折线 图 ; 实 线 为 拟 合 线 。 
图 7 年 均 蒸 散发 时 间 变 化 图 


Fig.7 Diagram of annual evapotranspiration changing 
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表 5 拟 合 回归 系数 


Tab.S Fit regression coefficient 


分 区 面积 /10* km? Bo B 
北 温带 区 (1) 18.73 -5517.40 2.80 
中 温带 区 (IT) 303.49 161.74 0.01 
南 温带 区 (IIT) 162.77 1645.51 -0.69 
北 亚热带 区 (IV) 56.76 -4028.99 2.27 
中 亚热带 区 (V) 120.64 -3788.96 2.14 
南亚 热带 区 (VD) 38.23 464.99 0.13 
北 热带 区 (VID) 7.60 310.23 0.27 
中 热带 区 (VIII) 0.91 -699.10 0.84 
高 原 气 候 区 (H ) 251.14 -268.73 0.27 
全 区 961.63 -1302.30 0.89 
注 :B 为 拟 合 回归 常数 项 ;8 为 拟 合 回归 一 次 项 系数 。 


表 6 中 国 不 同 分 区 T-S 分 析 系 数 和 M-K 分 析 Z 值 统计 
Tab.6 Statistics of T-S analytical index and M-K 


analytical Z value of each zone in China 


Ay T-S 分 析 M-K 分 析 ”双边 检验 

系数 Z 值 可 信和 度 /% 
北 温带 区 (1) 0.0055 1.54 87.64 
中 温带 区 (II) -0.0007 -0.03 2.40 
南 温带 区 (IITD) -0.0020 -0.40 31.82 
北 亚热带 区 (IV) 0.0048 0.51 39.00 
中 亚热带 区 (V) 0.0054 0.52 39.70 
南亚 热带 区 (VD) -0.0001 -0.02 1.60 
北 热带 区 (VII) 0.0003 0.01 0.80 
中 热带 区 (V1) 0.0015 0.11 8.76 
高 原 气 候 区 (H) 0.0003 -0.02 1.60 


注 :TS 分 析 系 数 为 Theil-Sen median 分 析 方 法 计算 的 趋势 斜率 ;M- 
开 分 析 Z 值 为 Mann-Kendall 分析 方法 计算 的 归 一 化 检验 统计 量 。 


4 结论 


<i 


(1) FP 22 EF By ais BUR. HH POE I8] ZR Pe A I 
大 ,主要 受 地 理气 象 要 素 影 响 。 各 分 区 蒸 散 发 中 中 
温带 区 (ID) . 南 温带 区 (IT) 、 北 亚热带 区 (IV) 和 高 原 
气候 区 (H) 空 间 变 异 系数 较 大 ,主要 原因 是 分 区 空 
间 上 的 异 质 性 较 高 ,其 他 分 区 变异 系数 较 低 且 空间 
地 理 环境 比较 稳定 。 

(2) 对 中 国 蔡 散发 的 时 间 特 征 分 析 显 示 ,在 相 
对 缺 水 干燥 的 北方 地 区 ,时 间 变 异 系数 一 般 较 大 ， 
而 南方 湿润 区 相对 系数 则 相对 偏 低 。 滨 海 城市 和 
盆地 边缘 山脉 附近 ,时 间 变 异 系数 较 大 , 受 不 规则 


0 500km 
Ls 


图 例 T-S 分 析 系 数 

i 通 量 站 。 0.00~0.22 
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审 图 号 : GS(2023)549 号 E 


注 :T-S 分 析 系 数 为 Theil-Sen median 分 析 方 法 计算 的 趋势 斜率 。 
图 8 2003 一 2020 年 中 国 工 S 分 析 系 数 空间 分 布 
Fig. 8 Spatial distribution of T-S analytical index 

in China from 2003 to 2020 
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注 :M-K 分 析 Z 值 为 Mann-Kendall 分 析 方 法 计算 
的 归 一 化 检验 统计 量 。 
图 9 2003 一 2020 年 中 国 M-K 分 析 Z 值 空间 分 布 
Fig.9 Spatial distribution of M-K analytical Z value in 
China from 2003 to 2020 


时 间 变 异 系数 差别 较 大 区 域 , 主 要 与 土地 覆盖 类 型 
BME TEAK. A KAR BUR TE TRIN Fe EB 
FE CS 2E, e DAE AR ROK ES To] ME HH BRUST 
间 并 不 完全 一 致 ,但 是 在 相关 邻 域 一 般 最 大 值 与 最 
小 值 出 现 的 时 间 较 为 相关 。 

(3) 通过 两 种 方法 得 出 北 温带 区 (1) 、 北 亚热带 


的 局 部 气候 条 件 影 响 较 大 。 同 时 ,各 分 区 内 也 存在 


KAV) 中 亚热带 区 (V) 、 北 热带 区 (VIL) .中 热带 区 
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图 10 2003 一 2020 年 中 国 赫 斯 特 指数 空间 分 布 
Fig. 10 Spatial distribution of Hurst index 
in China from 2003 to 2020 


(VII) 、 高 原 气 候 区 (H) 燕 散发 均 为 增长 趋势 , 南 温 
带 区 (II) 为 减 小 趋势 。 中 温带 区 (ID) 和 南亚 热带 区 
(VD) 两 种 分 析 方 法 所 得 的 变化 趋势 不 一 致 , 且 置 信 
度 水 平 较 低 ,主要 原因 可 能 是 两 区 蒸 散发 相对 稳 
AE ,波动 较 小 ,增长 或 减 小 趋势 不 明显 。 结 合 时 间 
变化 趋势 和 赫 斯 特 指数 对 未 来 中 国 蒸 散发 进行 预 
测 显 示 , 赫 斯 特 指数 除了 在 北 温带 区 (1) 大 于 0.5, 即 
未 来 的 变化 状况 与 过 去 趋势 一 致 ,但 赫 斯 特 指数 仿 
小 ,持续 性 较 小 。 其 他 区 域 赫 斯 特 均值 均 小 于 0.5， 
即 未 来 趋势 与 过 去 相反 ,但 值 与 0.5 接 近 ,持续 性 也 
不 强 。 本 次 研究 发 现 ,地 理气 候 因素 可 能 是 影响 大 
FUE DA ICY ETA RR 
(4) 人 研究 尺度 的 选择 对 于 问题 的 研究 息 息 相 
关 , 为 了 探究 影响 蒸 散发 在 时 间 空 间 上 的 不 均匀 性 
影响 因子 ,需要 通过 不 同 尺 度 ,不 同 分 析 数 据 进 
行 更 深层 次 的 研究 。 
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Temporal and spatial dynamic analysis of terrestrial evapotranspiration 
in China based on PML-V2 product 


WEI Tao, WANG Yunquan 
(School of Environmental Studies, China University of Geoscience(Wuhan), Wuhan 430074, Hubei, China) 


Abstract: Terrestrial evapotranspiration (ET) is an important process of land-atmosphere exchange and an impor- 
tant link in the global water migration and energy transfer system. Identifying the temporal and spatial dynamic 
characteristics of ET is of great significance in the study of regional water cycle and energy conversion. In this 
study, we compared the appropriateness of the three different ET products (GLEAM, MOD16, and PML-V2) in 
China using the flux data of nine tower stations in the country. We used the dataset of ET products between 2003 
and 2020 to analyze the temporal and spatial dynamic characteristics of evapotranspiration in China. The results 
showed the following: (1) The ET products of PML-V2 is the most suitable for China. (2) On the time scale, ET 
gradually increased during the research period. On the spatial scale, ET exhibited an increasing tendency during 
2003—2020 from northwest to southeast of China. (3) According to the Hurst index, the future ET of the whole 
nation is expected to show the opposite trend compared to the past. This means that ET increased in the past but 
is expected to decrease in the future. This study analyzed the spatiotemporal dynamic characteristics of ET in Chi- 
na, which could provide a reference for regional water resources utilization and optimal allocation. 


Key words: evapotranspiration; appropriateness; dynamic characteristic; China 


